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Computational Ion Optics Design Evaluations 


Shane P. Malone and George C. Soulas 
National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


Ion optics computational models are invaluable tools in the design of ion optics systems. 
In this study a new computational model developed by an outside vendor for use at the 
NASA Glenn Research Center (GRC) is presented. This computational model is a gun code 
that has been modified to model the plasma sheaths both upstream and downstream of the 
ion optics. The model handles multiple species (e.g. singly and doubly-charged ions) and 
includes a charge-exchange model to support erosion estimations. The model uses 
commercially developed solid design and meshing software to allow high flexibility in ion 
optics geometric configurations. The results from this computational model are applied to 
the NEXT project to investigate the effects of crossover impingement erosion seen during the 
2000-hour wear test. 




Nomenclature 



A 

= 

constant 

t 

= 

time 

B 

= 

constant 

T e 

= 

electron temperature, eV 


= 

accelerator grid aperture diameter 

T e d 

= 

downstream electron temperature, eV 

d s 

= 

screen grid aperture diameter 

T 

± eu 

= 

upstream electron temperature, eV 

e 

= 

elementary charge, 1.619xlO” 19 C 

v b 

= 

beam ion velocity, m/s 

E b 

= 

beam fast particle energy, eV 

Vid 

= 

downstream ion velocity, m/s 

ji 

= 

ion current density, A/m 2 

V a 

= 

accelerator grid voltage, V 

Jbeamlet 

= 

ion beamlet current density, A/m 2 

^net 

= 

net accelerating voltage, V 

h 

= 

ion beam current, A 

Vs 

= 

screen grid voltage, V 

m Xe 

= 

xenon atom or ion mass, 2.193xlCT 25 kg 

Xe 

= 

neutral xenon atom 

n b 

= 

beam fast particle density, nT 3 

Xe + 

= 

singly-charged xenon ion 

n e 

= 

electron number density, m -3 

Xe ++ 

= 

doubly-charged xenon ion 

Hi 

= 

ion number density, m -3 


= 

local plasma potential, V 

n id 

= 

downstream ion number density, nT 3 

<t>u 

= 

upstream plasma potential, V 

ftiu 

= 

upstream ion number density, nT 3 

<t>d 

= 

downstream plasma potential, V 

n n 

= 

neutral number density, nT 3 

ftp ’ 

= 

bulk plasma potential, V 

r 

= 

ion optics radial distance, mm 

®cex 

= 

charge exchange cross section, A 2 


I. Introduction 

Ion optics computational models are invaluable tools in the development of new ion engine optics and in the 
assessment of current ion optics designs. Ion optics models can be used as design tools to determine the best optics 
aperture size and grid spacing by assessing optics performance. A good ion optics model should be capable of 
predicting perveance of the optics, both direct and crossover impingement limits, beam divergence, and electron 
backstreaming limits. Ideally, the model should run relatively quickly and without the use of expensive mainframe 
computers. This allows the user to assess multiple iterations of a design in a short period of time. 

Ion optics models are also useful in assessing the service life and wearout mechanisms of a set of ion optics. By 
calculating the impingement of beam ions and charge exchange ions on the accelerator grid, ion optics models 
should be able to predict the erosion of ion optics over life. Performance modeling at different times in life with new 
wearout geometries can be used to assess the effects of erosion on perveance, beam divergence, and electron 
backstreaming. 
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Current and future ion thruster development programs will make extensive use of ion optics computational 
models to reduce design time and ion optics design iterations. This will allow optics to be designed and produced 
more effectively and quickly, at a lower cost to the development program. The NASA Evolutionary Xenon Thruster 
(NEXT) program has made use of ion optics models to determine the service life assessment of the ion optics, 1 ’ 2 and 
optics modeling was instrumental in the selection of Thick Accelerator Grid (TAG) optics geometry over the state- 
of-the-art NASA Solar electric propulsion Technology Applications Readiness (NSTAR) thruster ion optics 
geometry for the NEXT program. Additional ongoing work at NASA GRC in the development of high specific 
impulse, long life ion optics makes extensive use of ion optics modeling. 3 Other future thruster developments will 
likely require new ion optics designs and will also make extensive use of ion optics models to develop them. 

There are many ion optics models which have been developed in the last twenty years, some of which are 
described below. This is by no means a complete list of optics models, but is meant to illustrate the variety of 
approaches to the problem of ion optics modeling in the past and present. 

Ion optics models are most easily sorted into two-dimensional (2D) and three-dimensional (3D) models. Two- 
dimensional optics models are generally much faster to run than 3D models, owing to their relative simplicity and 
assumption of axisymmetry. Two-dimensional models have been developed and used at the University of 
Michigan, 2 University of Tokyo, 4 and the Keldysh Research Center. 5 Unfortunately, the assumption of axisymmetry 
which allows the 2D models to be run much faster than 3D models introduces some inherent limitations in their 
capabilities, namely that they cannot predict over-focusing well, and that they cannot determine impingement 
patterns on downstream surface of the accelerator grid in a two-grid optics system. However, for situations where 
barrel erosion is the dominant wearout mechanism, such as three-grid optics systems and some two-grid optics 
systems, 2D codes have applications in service life assessment. 

Three-dimensional optics models are generally more useful than 2D models for optics design and service life 
assessment, with the drawback that many models cannot be run quickly on personal computers without some 
simplifying assumptions. This limits the amount of design iterations that can be evaluated. Three-dimensional 
models have been developed by University of Tennessee Space Institute, 6 Colorado State University, 7 University of 
Tokyo, 8 and jointly by the Virginia Polytechnic Institute and the Jet Propulsion Laboratory. 9 

In this paper, a commercially developed, fully three-dimensional ion optics computational model with 
tremendous flexibility is presented. This model uses commercial, off-the-shelf computer aided design software and a 
similarly available commercial meshing package as the front end, allowing modeling of single and multiple 
apertures, varied aperture shapes and grid spacing, and the commonly seen cusp geometry along grid aperture 
barrels, as well as misaligned apertures. This geometric flexibility allows investigations of new geometrical 
concepts. The model is also capable of modeling multiple charge state ions and multiple generations of charge 
exchange (CEX) ions and neutrals. Through the use of reflective boundary conditions, the model can simulate 
adjacent apertures rather than modeling them directly, and can also simulate the effects of edge apertures and 
“missing” holes along the edge where the hole pattern ends. This optics model is currently being used to evaluate the 
TAG optics on the NEXT thruster. Preliminary results from this evaluation are presented. 

II. Ion Optics Computational Model Description 

The ion optics computational model developed for NASA GRC is a new three-dimensional electrostatic 
equilibrium particle-in-cell (PIC) code that has been designed to address the recent beam optics modeling and 
simulation requirements for vacuum electron devices, ion sources, and charged-particle transport. The model 
combines modern finite-element techniques with improved physics models. The code employs a conformal mesh, 
including both structured and unstructured mesh architectures for meshing flexibility, along with a new method for 
accurate, efficient particle tracking. New particle emission/creation models include space-charge-limited emission, 
charge exchange, thermionic beam representation, and secondary electron emission. One of the key features of the 
code is its ability to model fine-scale features in a large volume. This makes it valuable for performing sensitivity 
studies and generating manufacturing tolerance specifications. 

Although the code is still in development, it has been specifically designed to model ion thruster optics. 
Currently, much emphasis is being placed on validating the code against experimental data and other models. This is 
especially necessary for ion thruster modeling, where detailed knowledge of the physical parameters is difficult to 
ascertain. 

For ion thruster modeling, the need for three-dimensional simulation of volumetric ion sources, ion acceleration, 
and optics geometry, with the ability to model charge exchange of the ion beam with a background neutral gas are 
important features to be captured. The two pieces of physics that stand out as significant are modeling of the 
upstream and downstream plasma sources and modeling of charge exchange. 
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The numerical solution employed by the code combines finite element analysis with particle-in-cell (PIC) 
methods. The basic physics model is based on the equilibrium steady-state application of the electrostatic particle- 
in-cell (PIC) approximation employing a conformal computational mesh. The foundation of the model stems from 
the same basic electron gun code model introduced in codes such as EGUN. 10 Here, Poisson's equation is used to 
self-consistently include the effects of space charge on the fields, and the relativistic Lorentz equation is used to 
integrate the particle trajectories through those fields. 

Gun models have been around since the 1970s, and have been a mainstay in the vacuum electronics and 
accelerator communities for beam formation modeling. As shown in Figure 1, in the gun model the fields are held 
frozen while the particles are tracked through the system. As the particles are tracked, they lay down charge on the 
mesh, which is subsequently used to update the fields. This series of steps is repeated until the beam tracking and 
fields stop evolving, and results in the self-consistent steady-state equilibrium condition. The benefit of this method 
over time-domain models is that there are no time-related transients that have to be flushed out of the simulation, 
and the expensive field solution step is not done every step. These simulations can be orders of magnitude faster 
than time-domain solutions. 

For ion thruster modeling, several modifications and additions to the basic electron gun code model needed to be 
made, including: 

• inclusion of ions and neutrals, 

• a model to support the ion plasma source, and 

• addition of a charge exchange model. 

These new features were not trivial to add to the gun model, but fortunately the gun model supports the architecture 
needed for their implementation. 

In the new code, the method for modeling the plasma sheath that forms in ion sources assumes that the electron 
distribution function is a Maxwellian function of electrostatic potential over electron temperature. The plasma 
model feature supports the development of a sheath in both the upstream ion source and downstream plasma. 
Separate upstream and downstream models are included that supports thruster modeling. The electrons are included 
numerically as Boltzmann electrons with a density calculated per Equation 1. 


n 


e 



( 1 ) 


The source faces of the mesh can support multiple species per source face so that Xe + and Xe ++ can be injected 
simultaneously. The model is not constrained by imposing a Dirichlet boundary condition that dictates the potential 
at the plasma boundaries. As a result this potential is calculated self-consistently. However, the user can impose a 
Dirichlet boundary condition if desired. 

Charge exchange is the process by which a “fast” charged particle streaming through a neutral background gas 
collects an electron from a “slow” neutral particle, resulting in a fast neutral and a slow charged particle. An 
efficient method for capturing this is essential. The model supports multiple charge exchange reactions, including 
the following: 


Xe+ fas, + Xe slow Xe fast 


+ Xe 


slow 


Xe 


++ 

fast 


+ Xe slow 


- Xe fast 


+ Xe 


++ 

slow 


Xe fast + Xe slcw 


X f a ,+ X e 


slow 


The charge exchange collision rate follows Equation 2. A semi-empirical model by Miller, et al 11 was used for 
the charge exchange algorithm, shown in Equation 3. The CEX algorithm implemented handles multiple 
generations of charge exchange. The self-consistent space charge is included in the Poisson solution, so that the 
effects of source and CEX products are captured. The subsequent fast neutral particles are tracked for possible 
collisions with the grids. The model currently performs one “reaction” per cell, where the source particles are 
binned and multiple launches can be made per “reaction”. Here, the source ions are velocity-binned (bins chosen by 
user), and each launch creates one launch pair per bin. This allows the code to capture the effect of a distribution, 
while allowing the total particle count to stay in check. 


dn b 

dt 


( 2 ) 
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Ocex = A- B lo g(E b ) 


( 3 ) 


III. Model Application to NEXT 

A. Problem 

During the NEXT 2000-hour wear test, accelerator apertures in the outer 47 mm of the grids exhibited an off- 
center, hexagonal star-shaped erosion pattern. 12 The eroded points of the stars were chamfered through the electrode 
thickness and extent of the erosion pattern generally increased with increasing ion optics’ radius (see Figure 2). It 
was determined that this erosion likely occurred early in the wear test and at lower beam current operation (1.20 A 
and 2.70 A). Erosion was most likely due to a combination of effects, including: 

1. over-focusing of beamlet ions due to low plasma density near the grid edges; 

2. significantly smaller-than-nominal accelerator aperture diameters in this region; and 

3. misaligned screen and accelerator apertures due to an under-compensated hole pattern. 

Over-focusing of ion beamlets generally occurs in areas where there is low plasma density, such as near the 
edge apertures and when the thruster is operated at a lower beam current. The normalized perveance per hole drops 
to values much less than one, and the plasma sheath migrates further upstream. This highly distorted sheath tends to 
draw ions across the beam centerline, where, in extreme cases, they impinge on the accelerator grid aperture barrel. 
This effect is generally referred to as crossover impingement. 

Due to manufacturing difficulties associated with etching the TAG optics on the NEXT engine, hole sizes near 
the edge were more than 10% smaller than the nominal hole size. Figure 3 and Figure 4 show the accelerator and 
screen grid aperture sizes as a function of radius. Smaller holes in the accelerator grid mean that it is easier for over- 
focused ions to strike the barrel of the aperture. 

Under-compensation of the hole pattern means that the screen and accelerator grid apertures are not co-linear 
but are misaligned such that the outer beamlets are deflected away from the thrust axis. A misaligned set of 
apertures would be more susceptible to grid erosion due to electric fields which pull the ion beamlet away from the 
geometric centerline of the holes. 

B. Optics Model Assessment 

The ion optics computational model is being used to assess the over-focusing of the beamlets in the low density 
regions of the NEXT ion optics. Preliminary evaluation of over- focusing is being performed for engine operation at 
3.52 A beam current and 1.20 A beam current, at radial locations of 153 mm (where erosion patterns were first 
noticed) and 199 mm (the last edge aperture). Current density values and aperture sizes are directly from 
measurements taken during the NEXT 2000-hour wear test and subsequent examination. Aperture sizes were 
determined by averaging the positive and negative aperture sizes shown in Figure 3 at each radial location. 
Thruster current density profiles for 3.52 A and 1.20 A cases are shown in Figure 5. A nominal grid was used in 
these calculations. Model inputs are shown in Table 1. 

Screen and accelerator voltages were those used during the 2000-hour test. Plasma parameters were determined 
based on plasma diagnostics during the test. Upstream plasma density was determined by assuming that the 
upstream plasma ion current was the Bohm current, and calculated from Equation 4. 



The downstream plasma density was estimated using Equation 5, where ion velocity downstream is determined 
from conservation of energy, Equation 6. 


( 5 ) 


( 6 ) 
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The effects of CEX ions and neutrals were not considered for this initial analysis, but will be investigated at a 
later time. 

C. Geometric Configurations 

Several model mesh geometries were used to attempt to determine the effects of over-focusing and crossover 
impingement. These included co-linear screen and accelerator grid apertures using both smooth barrels and the cusp 
geometry shown in Figure 6, as well as apertures which were misaligned by up to 4% of the nominal screen grid 
aperture diameter, as shown in Figure 7. 

In most cases, a hexagonal, reflective boundary condition (see Figure 8) was used to determine the proper ion 
impingement patterns on the accelerator grid aperture barrel, but a cylindrical boundary condition was used for cusp 
cases because the model was easier to create. Incorporation of the cusp geometry with the hexagonal reflective 
boundary condition will be implemented in the future. 

Because areas of the ion thruster with low plasma density and low ion current density are being investigated, a 
long upstream area is required in order to allow the upstream plasma potential to reach equilibrium and the sheath to 
form properly. The upstream area length is 3.2 times longer than the nominal screen grid aperture diameter. This 
provides adequate upstream distance for the extended sheath and pre-sheath to set up and reach equilibrium at the 
plasma potential. In order to reduce computational time, mesh elements at the upstream boundary to the model are 
quite long, but get smaller as they approach the screen grid. This allows the model to focus more on the strong 
electric field region of the optics. The downstream area is similarly configured, with mesh elements growing 
longer as they get further from the downstream surface of the accelerator grid. An example of the mesh is shown in 
Figure 9. 

There are 1947 source faces at the upstream mesh boundary. From each source face, 3 emission macroparticles 
are injected. Each macroparticle represents a number of ions such that the total current density required for the case 
is maintained. 

D. Preliminary Results 

Figure 10 and Figure 11 show the equipotential lines and ion trajectories for the r = 199 mm cases. Normalized 
perveance per hole was calculated per Reference 3, giving values of 0.006 and 0.021 for the J b = 1.20 A and 3.52 A 
cases, respectively. The J b = 1.20 A case shows clear evidence of crossover impingement along the accelerator grid 
aperture barrel by over-focused beam ions. However, the J b = 3.52 A case shows that there is either no impingement 
or very little impingement, a result that is surprising given the extent of the erosion shown in Figure 2. Further 
refinement of the model geometry and input parameters will be necessary to determine whether crossover 
impingement does actually occur for this operating condition. 

For the J b = 1 .20 A case, impingement current density was plotted on the 3D surface of the accelerator aperture 
barrel to determine the qualitative erosion pattern (see Figure 12). A future improvement in the model will include 
the ability to show impingement current for ions of different energies in order to calculate the erosion rate directly. 
The most current is deposited near the downstream surface of the accelerator aperture barrel, and in six directions 
indicating the points of the hexagonal boundary. The points of the hexagonal boundary correspond to the webbing 
directions as shown in Figure 8. As shown in Figure 2, for apertures which are fully surrounded by holes (r = 184 
and 195 mm pictures), there is preferential erosion along the webbing direction. This is most likely due to upstream 
sheath distortion caused by the proximity of adjacent holes. Figure 2 shows that at r = 199 mm, where the hole 
pattern ends, a different erosion pattern emerges, most likely due to the “missing” holes beyond the hole pattern. 
This effect will be investigated further in the future. 

Figure 13 and Figure 14 show the equipotential lines and ion trajectories for the r = 153 mm cases. Normalized 
perveance per hole was calculated per Reference 3, giving values of 0.055 and 0.161 for the J b = 1.20 A and 3.52 A 
cases, respectively. Neither the J b = 1.20 A nor the J b = 3.52 A cases show any evidence of crossover impingement. 
This most likely indicates that the erosion seen in the NEXT 2000-hour test at this aperture locations was not caused 
primarily by the low upstream plasma density as in the r = 199 mm case. To investigate other geometric reasons for 
the erosion, the effects of cusps and misaligned apertures were investigated. The model was run for these geometries 
using the same conditions shown in Table 1 for the J b = 1.20 A case. 

Figure 15 and Figure 16 show the equipotential lines and particle trajectories for the cusp and misaligned 
aperture cases, respectively. There is still no evidence of impingement current on the accelerator aperture barrel for 
the cusp geometry. However, the electric field lines are distorted differently due to the change in geometry, an effect 
which could cause changes in the ion trajectory. Future calculations will combine the hexagonal boundary with the 
cusped aperture to provide the most realistic look at the optics geometry. 
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The misaligned apertures create electric fields which tend to slew the ion trajectories away from the direction of 
misalignment (in the case of the NEXT 2000-hour test, away from the thrust axis). This, combined with the low 
upstream plasma density, creates an over-focused beam which impacts the aperture barrel preferentially as shown in 
Figure 16. This preferential impingement is likely the cause of the slight erosion shown in Figure 2 for the r = 153 
mm case. 


IV. Future work 

The ion optics computational model being developed for NASA-GRC will be used in the design of ion optics for 
future ion thrusters. Further investigation of the NEXT 2000-hour wear test results will include: 

• Incorporation of the cusp geometry with the hexagonal boundary condition 

• Investigation of varied grid gap dimension on crossover impingement 

• Modeling of edge apertures to simulate “missing” apertures at the edge of the hole pattern 

• Sensitivity analysis of input parameters such as T eup and (f up on crossover impingement 

• Inclusion of CEX effects, including pit-and-groove erosion analysis 

• Incorporation of ion energy with impingement current density to determine erosion 

In order to more fully validate the results of the ion optics model, it will be compared to the results of the 
NSTAR LDT and ELT, as well as further data from the NEXT 2000-hour wear test. The model will also be 
compared to the results of other established computational models to validate certain characteristics. 

V. Conclusions 

The ion optics computational model being developed for NASA-GRC shows great flexibility in geometry due to 
its use of commercial computer aided design and meshing software as the front end. The code itself is capable of 
running on a personal computer, with detailed runs taking only a few hours to finish. The model was used to 
investigate crossover impingement erosion seen in the NEXT 2000-hour wear test. Preliminary comparison of 
model results to the wear test results demonstrate the model’s flexibility in modeling various aperture geometries, 
but indicate that improvements to the model and to understanding of the input conditions must be made to fully 
predict over-focusing and crossover impingement. Further improvements in the model will enhance this capability. 
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Table 1. Model inputs for NEXT over-focusing. 



J b = 

r = 153 mm 

1.20 A 

r = 199 mm 

Jb = 

r = 153 mm 

3.52 A 

r = 199 mm 

da 

-6.7% from 

-15.6% from 

-6.7% from 

-15.6% from 


nominal 

nominal 

nominal 

nominal 

d s 

-2.7% from 

-2.7% from 

-2.7% from 

-2.7% from 


nominal 

nominal 

nominal 

nominal 

V a 

-260 V 

-260 V 

-220 V 

-220 V 

v s 

1765 V 

1765 V 

1766 V 

1766 V 

</>u 

1796 V 

1796 V 

1796 V 

1796 V 

</>d 

6 V 

6 V 

6 V 

6 V 

T eu 

6 eV 

6 eV 

6 eV 

6 eV 

Ted 

1 eV 

1 eV 

1 eV 

1 eV 

n iu 

3.30xl0 16 m“ 3 

3.38xl0 15 m -3 

9.89xl0 16 m“ 3 

1.29xl0 16 m -3 

n id 

1.37xl0 15 m“ 3 

1.41xl0 14 m -3 

4.12xl0 15 m -3 

5.37xl0 14 m -3 

Jbeamlet 

11.18 A/m 2 

1.14 A/m 2 

33.54 A/m 2 

4.37 A/m 2 
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Figure 1. Basic gun code algorithm 



r = 195 mm r = 199 mm 


Figure 2. Post-test photomicrographs of NEXT downstream accelerator apertures at radial locations of 153, 
184, 195 and 199 mm. Aperture erosion is less obvious at 153 mm but appears in the 4 o’clock area of the 
aperture. Aperture erosion at 184 and 195 mm shows preferential erosion along the webbing direction. 
Aperture erosion at 199 mm is distorted due to the absence of apertures to the right of the aperture shown. 
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Accelerator Grid Aperature Diameter (% of Nominal) 



Radial Location (mm) 

Figure 3. Pre-test accelerator aperture diameters as a function of radius. Diameters were measured with pin 
gages that were in increments of 2.2% of the nominal accelerator aperture diameter. 
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Figure 4. Pre-test screen aperture diameters as a function of radius. Diameters were measured with pin gages 
that were in increments of 1.3% of the nominal screen aperture diameter. 
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Figure 5. Beamlet current density as a function of radius for 1.20 A and 3.52 A beam current conditions. 

Line indicates location where aperture erosion begins. 


NAS A/TM— 2004-2 1 3206 


10 



D* 

Figure 6. Aperture cusps were modeled. Aperture diameter at the cusp tips is the measured value in Figures 
2 and 3. Cusp height is 20% of the nominal screen grid thickness and 33% of the nominal accelerator grid 

thickness. 



z 

Px 

Figure 7. Misaligned apertures were modeled. Accelerator aperture centerline is offset from screen aperture 
centerline by 4% of nominal screen grid aperture diameter. 


Boundary 



Modeled 

Aperatures 


Figure 8. The ion optics model uses a full single aperture and a hex shaped reflective boundary to simulate 

adjacent apertures. 
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Figure 9. A cross-section of the hexagonal mesh geometry is shown in (a). A lengthwise cross-section of the 

mesh geometry is shown in (b). 
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Figure 10. Equipotential lines and particle trajectories for r = 199 mm, J b = 1.20 A case. The sheath is 
distorted well upstream of the screen grid, and crossover impingement is clearly indicated. Boundaries 
shown are from point to point on the hex boundary condition. 
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Figure 11. Equipotential lines and particle trajectories for r = 199 mm, J b = 3.52 A case, 
distorted upstream of the screen grid, but crossover impingement is not clearly indicated. 

are from point to point on the hex boundary condition. 
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Figure 12. Impingement current on the accelerator grid barrel and downstream face due to crossover 
impingement for r = 199 mm, J b = 1.20 A case. Notice that the highest impingement current is deposited in 
the directions of the hexagon points, which correspond to the webbing direction of the inter-aperture space. 
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Figure 13. Equipotential lines and particle trajectories for r = 153 mm, J b = 1.20 A case. Notice that sheath is 
still established relatively near the screen grid in this case. No crossover impingement is indicated. 
Boundaries shown are from point to point on the hex boundary condition. 
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Figure 14. Equipotential lines and particle trajectories for r = 153 mm, J b = 3.52 A case. The sheath is 
again established close to the screen grid. No crossover impingement is indicated. The apparent asymmetry 
in the ion trajectories is most likely due to insufficient run cycles, and not an actual physical effect. 
Boundaries shown are from point to point on the hex boundary condition. 
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Figure 15. Equipotential lines and particle trajectories for r = 153 mm, J b = 1.20 A, cusp geometry case. The 
cusps distort the equipotential lines inside the aperture barrels, but no crossover impingement is indicated. 

The boundary condition is cylindrical (axisymmetric). 
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Figure 16. Equipotential lines and particle trajectories for r = 153 mm, J b = 1.20 A, misaligned aperture case. 
The sheath is still distorted upstream of the screen grid, and crossover impingement is clearly indicated on 
the side of the aperture barrel opposite of the misalignment direction. Boundaries shown are from point to 

point on the hex boundary condition. 
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